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3 ■ ABSTRACT 

\ 

^ : 

^ ■ The ever-expanding depth and quahty of photometric and spectroscopic ob- 

O ■ servations of stellar populations increase the need for theoretical models in regions 

■ of age-composition parameter space that are largely unexplored at present. Stel- 
^ ■ lar evolution models that employ the most advanced physics and cover a wide 

■ range of compositions are needed to extract the most information from current 
observations of both resolved and unresolved stellar populations. The Dart- 
mouth Stellar Evolution Database is a collection of stellar evolution tracks and 
isochrones that spans a range of [Fe/H] from -2.5 to +0.5, [a/Fe] from -0.2 to 
+0.8 (for [Fe/H] <0) or +0.2 (for [Fe/H] >0), and initial He mass fractions from 
Y=0.245 to 0.40. Stellar evolution tracks were computed for masses between 0.1 
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and 4 allowing isochrones to be generated for ages as young as 250 Myr. For 
the range in masses where the core He flash occurs, separate He-burning tracks 
were computed starting from the zero age horizontal branch. The tracks and 
isochrones have been transformed to the observational plane in a variety of pho- 
tometric systems including standard UBV(RI)c, Stromgren uvby, the Sloan Dig- 
ital Sky Survey ugriz, the Two Micron All Sky Survey JHKg, and Hubble Space 
Telescope ACS-WFC and WFPC2. The Dartmouth Stellar Evolution Database 
is accessible through a websit^ where all tracks, isochrones, and additional flies 
can be downloaded. 

Subject headings: globular clusters: general — open clusters: general — stars: 
evolution 



1. Introduction 

Recent photometr ic surveys such as the ACS Survey of Galactic Globular Clusters 
( ISarajedini et al.ll2007l ) have substantially increased the quality and volume of data available 
for the Galactic globular cluster system. Despite the fact that this survey focused on metal 
poor stellar populations it was found that there was a need for stellar evolution models with 
[Fe/H] >0 (jSiegel et al.l 120071 ). Many other photometric surveys such as the Two Micron 
All Sky Survejo (2MASS) and the Sloan Digital Sky Survej|§ (SDSS) and extensions such 
as SEGUeQ underscore the need for stellar evolution libraries that span a wide range of 
compositions and photometric systems. 



Dotter et al.l (120071 ) introduced a set of stellar evolution tracks, isochrones, and computer 



programs to accompany the ACS Survey of Galactic Globular Clusters. The set includes 
tracks with masses between 0.1 and 1.8 M©, isochrones with ages ranging from 2 to 15 Gyr, 
[Fe/H] from -2.5 to with a variety of options in [a/Fe] and initial He content {Yinu)- In 
this paper, that set of models is expanded to include: (i) tracks and isochrones with [Fe/H] = 
-1-0.15, -1-0.3, and -1-0.5; (ii) higher mass tracks and isochrones with ages as low as 250 Myr for 
all compositions; (iii) a wider selection of photometric systems; and (iv) a web interface to 



^http: / /stellar. dartmouth.edu/~models/ 
^ lit tp : / / pegasus . phast . umass . edu/ 
'^http://www. sdss.org/ 
^http://segue. uchicago.edu/ 
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three of the computer programs so that visitors can easily generate specific models without 
downloading the entire database and running the programs locally. 

Collectively, these tracks, isochrones, and supporting computer programs constitute the 
Dartmouth Stellar Evolution Database. The stellar evolution code used to generate these 
models is reviewed in §2. Details of the compositions and model grids are outlined in §3 
and the color-Tefr transformations in §4. In §5, the isochrones are compared to those of 
comparable isochrone libraries and photometry of Galactic open clusters and one globular 
cluster. The results are summarized and discussed in 56. 



2. The Dartmouth Stellar Evolution Program 

The Dartmou t h Stel lar E volution Prog r am ( DSEP) has been described previously by 
Biork fc Chaboveil J2006h and botter et all j2007h . The following is a brief description of 
the physics employed by DSEP. 

Partially inhibited atomic dif fusion and gravit a tiona l settling of both He (Y) and metals 
(Z) are included as described by Chaboyer et al. (j2001 ). Diffusion is calculated based on 
the formalism of iThoul et al.l (Il994j ) but inhibited in the outer 0.01 Mq where the diffusion 
coefficients are set to zero in the outermost 0.005 Mq and ramped from zero to full diffusion in 
the inner 0.05 Mq. This m ethod produces re sults that are consistent with recent observations 
of diffusion in NGC 6397 korn et al.ll2006h . 



DSEP uses high temperat ure opacities from OP AL (llglesias fc Rogers! Il996l ) and low 
temperature opacities based on [Ferguson et al.l (120051 ). The OPAL tables were obtained for 
each of the different [a/Fe] mixtures from the OPAL web interface § and the low temper- 
ature tables were calculated by one of us (Ferguson). The low temperature opacity tables 
have been made available via the World Wide Weqj. Conductive opacities are taken from 



Hubbard fc Lampd (jl969|). 



The basic equ ation of state (EOS) in DSEP is a general ideal gas EOS with the Debye- 
Hiickel correction (jChaboyer fc Kimlll995l ) and is used for t racks with M > O.SM©. For lower 
mass tracks (M< O.SMq), the FreeEOS equation of state (Ilrwinll2004l ) is used in the E0S4 
configuration. 

DSEP treats convection with the standard mixing length theory and a solar-calibrated 



'"^http: / /www-phys. llnl.gov/Research/OPAL/ 
®littp:/ /webs. Wichita. edu/physics/opacity/ 
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mixing length, a/vff,= l-938. The mode l s inclu de convective core overshoot according to the 
scheme developed by iDemarque et al.l (120041 ) . The underlying idea is that the amount of 
overshoot grows with the size of the convective core. The amount of core overshoot is 
parametrized as a multiple of the pressure scale height (Ap) and is a function of stellar mass 
and composition. The minimum stellar mass receives 0.05 Ap of overshooting, models with 
masses greater than the minimum by 0.1 Mq receive 0.1 Ap of overshooting, and models with 
masses greater than or equal to the minimum plus 0.2 Mq receive 0.2 Ap of overshooting. For 
details on the minimum mass as a function of composition see Table [H Convective envelope 
undershooting is not considered in these models. 

The su rface boundary conditio n (BC) was derived from a grid of PHOENIX model at- 
mospheres (IHauschildt et al.l ll999al Jbl) covering Teff=2,000 to 10,000 K and log g=-0.5 to 
5.5. The solar abundance scale used in these calculations was that of Grevesse fc Sauval 



(Il998l . hereafter GS98). PHOENIX model atmosphere grids were computed for the full range 
[Fe/H]=-2.5 to +0.5 and [a/Fe]=-0.2 to +0.8 covered by the stellar evolution calculations. 
A small number of model atmospheres with enhanced He were also computed to gauge the 
importance of enhanced He to the physical structure and the spectrum. The gas pressure 
(and thus the surface BC) experiences a small change with enhanced He and this is accounted 
for in the surface BC. 



Castelli &: Kuruczl (120031 ) atmospheres were used for T^s > 10,000 K; ICastelli fc Kurucz 



models were computed only for [Q;/Fe]=0 and +0.4. Stellar evolution models with 
[«/Fe] < +0.2 made use of surface BC derived from [a/Fe] = model atmospheres for 
Tcff > 10,000 K. For stellar evolution models with [a/Fe] > +0.4 the [a/Fe] = +0.4 bound- 
ary conditions we r e used for Tcfj > 10,000 K. The transition between the PHOENIX and 
Castelli fc Kuruca (120031 ) were always smooth regardless of the composition. The smooth 
transition is due to the stellar evolution models being relatively insensitive to the surface BC 
at higher Tcfr coupled with the reduced sensitivity of the model atmospheres to abundance 
variations in the transitional T^s- 

The surface BC was fit at the point in the atmospheric structure where T=Tefj for a 
given model. The gas pres sure was extracted at T= Tpff and added to the radiation pressure 
to satisfy the surface BC. IVandenberg et al.l (120071 ) found that, at least in the vicinity of 1 
Mq, the choice of model atmosphere boundary condition fitting point (either at T=Tcff or 
at an optical depth of r=100) has little infiuence on the resulting evolution. 

The PHOENIX model atmosphere code employs its own set of physics that is not identical 
to the physi cs used in DSEP , thoug h there is some overlap. In particular, the low temperature 
opacities of iFerguson et al.l (120051 ) are calculated with PHOENIX and this implies a level of 
consistency between the two codes. PHOENIX employs a standard mixing length theory of 
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convection with mixing length umlt = 2.0 (whereas DSEP uses umlt = 1-938). Perhaps 
the largest difference between PHOENIX and DSEP in terms of stellar structure is the EOS: 
PHOENIX uses an ideal EOS that includes the contributions of a large number of elements and 
molecules while FreeEOS (used by DSEP) includes several non-ideal effects but accounts for 
fewer elements and only molecules involving combinations of H and He. 



The nuclear reaction rates are taken fr om lAdelberger et al.l (119981 ) for most reactions. 
The exceptions are the ^^N (p,7)^^0 react i on (|lmbriani et al.ll2004j ). the triple-a r eaction from 



the N ACRE compilation (lAngulo et al.lll999f ) ; and the ; ^ ^C(q;,7 )^^0 reaction (IKunz et al. 



2OO2I ). Neutrino cooling rates were taken from lHaft et al.l (119941 ). 



3. Model Grids and Supporting Programs 

Stellar evolution models were computed using one of six different heavy element mix- 
tures with the distinguishing characteristic of each mixture being the level of a-element 
enhancement. Opacities were generated for each of [a/Fe]= -0.2, 0, 0.2, 0.4, 0.6, and 0.8. 
The heavy element mass fractions of these mixtures are given in Table [2l The mixtures 
were calculated starting with the solar mixture of GS98 by adding the value of [a/Fe] to the 
log(N-|-12) abundance of each of the following a-elements: O, Ne, Mg, Si, S, Ca, and Ti; the 
abundance of Ar was not altered from its GS98 value in any of the mixtures. 

Overall, models were computed for 67 different X,Z pairs; the initial compositions are 
listed in Table [31 For [Fe/H] > the maximum level of a-enhancement was chosen to 
be [a/Fe]=+0.2 in order to keep the total Z value below 0.1 since the opacity tables do 
not extend above Z=0.1. The model atmosphere grid included only provided scaled-solar 
([a;/Fe]=0) models at [Fe/H]=+0.1 5. The He content of the models begins with Y=0.245 
at [Fe/H] = -2.5, consistent with ISpergel et al. J2003h . and a He enrichment relation of 



AY/AZ=1.54. For [Fe/H] <0, additional models were computed with Y=0.33 and 0.40 (see 
Table ED . 

Stellar evolution tracks were computed for stellar masses ranging from 0.1 to 4 in 
increments of 0.05 Mq (between 0.1 and 1.8 M©), 0.1 M© (between 1.8 and 3 Mq), and 0.2 
Mq (between 3 and 4 Mq). The tracks cover the evolution of the models from the fully 
convective pre-main sequence (pre-MS) to either the core He flash for M < 2 Mq or either 
100 Gyr or the top of the white dwarf cooling curve (whichever came first) for M < 0.5 M©. 
For models that experience the core He flash, separate He burning tracks were evolved from 
the Zero Age Horizontal Branch (ZAHB) using the surface compositions and He core masses 
listed in Tables H] and O respectively. For models that either transition smoothly to core 
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He burning (M > 2 Mq) or were evolved from the ZAHB, the He-burning evolution was 
followed until the onset of thermal pulsations on the asymptotic giant branch (AGB). 

Isochrones were generated from these tracks with ages from 250 Myr to 15 Gyr in 
increments of: 50 Myr between 250 Myr and 1 Gyr, 250 Myr between 1 and 5 Gyr, and 
500 Myr between 5 and 15 Gyr. The isochrones extend from the lower MS to the tip of 
the red giant branch (RGB) or, for ages where the tracks permit it (MS mass above ~2 
Mq), at or near the onset of the thermal pulsations on the AGB 0. The low mass limit of 
the isochrones depends on the age: the isochrones only include masses for which the models 
have reached the MS. The lowest mass stars considered in this paper (0.1 M©) may require 
~1 Gyr before reaching the MS and therefore will not be present in isochrones for younger 
ages. For example, the [Fe/H]=0 isochrones have the following minimum masses (Mq): 0.43 
at 250 Myr, 0.32 at 500 Myr, 0.22 at 750 Myr, and 0.11 at 1 Gyr. Evolutionary tracks are 
also provided and these include the pre-MS evolution for all masses. 

All tracks and isochrones were transformed to the color-magnitude diagram (CMD) 
using each of the color-Tefj transformations described in §4. The track files consist of lines 
indicating age, log Tefr, log g, log L/Lq, followed by the absolute magnitude in each of the 
filters for a given photometric system. A theoretical track file is also available, this file 
contains age, log T^s, log g, log L/L©, log R/R©, core He mass fraction (Ycore), core Z mass 
fraction (Zcore); the surface Z/X ratio, H-burning luminosity (Lq), He-burning luminosity 
(Lq), He core mass, and C/0 core mass. The isochrones list stellar mass, log T^s, log g, log 
L/L0, and absolute magnitude in each of the filters for a given photometric system. 



The collection of Fortran 77 programs introduced by iDotter et al.l (120071 ) have been 



updated to deal with the expanded compositions and photometric systems presented i n this 



paper but are otherwise largely unchanged. See the appendices of iDotter et al.l (120071 ) for a 
detailed description of these programs and their output. 



4. Color-Teff Transformations 

Synthetic color-Teg transformations were derived from the PHOENIX model atmosphere 
grid (also used to generate surface boundary conditions in the stellar evolution calculations) 



^Due to the mass spacing of the stellar evolution tracks and the interpolation used in isochrone generation, 
the isochrones that include the AGB do not always terminate at the 'tip' of the AGB. In this case, it is 
recommended to examine the stellar evolution tracks if the brightness of the stars at the onset of the TP- 
AGB stage is of interest. This applies only to iscohrones with ages below 1 Gyr. Isochrones with ages above 
1 Gyr that include only the RGB provide reasonably accurate estimates of the RGB tip as a function of age. 



- 7- 



for Tcff < 10,000 K and for all [Fe/H] an d [a/Fe] combinations. The PHOENIX grid was 
supplemented by ICastelli fc Kuruca (120031) mod els for higher Teff. In order to insure a 
smooth transition, the ICastelli fc Kuruczl (120031 ) colors were adjusted to hne up with the 
PHOENIX colors at 10,000 K and the two sets of tables were ramped between PHOENIX and 
CasteUi fc Kuruczl J2003h in the temperature range of 9,000 K to 10,000 K. 



As mentioned in §2, the PHOENIX model atmospheres used to generate synthetic color- 
Tefj transformations were computed for each of the [Fe/H], [a/Fe] pairs but not for all of the 
enhanced Y values. The influence of enhanced He on broad-band synt hetic photoraetry is 
negligible for the tem peratures cons i dered here, as discussed briefly by iDotter et al.l (120071 ) 
and in more detail by iGirardi et al.l (120071 ). 



The ICastelli &: Kuruczl (120031 ) model atmospheres were only computed for [a/Fe]=0 and 
+0.4. In constructing the color transformations the same approach was used as described 
in §2 for the surface BC. Fo r evolu tionary tracks and isochrones with [a/Fe] < +0.2, colors 
based on ICastelli fc Kuruczl (120031 ) synthetic spectra wi th [g/Fe] = were applie d and for 



tracks and isochrones with [a/Fe] > +0.4, colors based on lCastelli fc Kuruczl (120031 ) synthetic 
spectra with [a/Fe] = +0.4 were applied. While this is not an ideal situation, the influence 
of a-enhancement, or indeed of total Z, on broad-band synthetic colors is small for Tcs ^ 
5,000 K. The only exceptions are in the bluest bands (equivalent to U) and at Tefr > 35,000 
K. Bands equivalent to U are the most uncertain in terms of the synthetic photometry and 
very few of the evolutionary tracks (only the lowest mass, metal poor He-burning tracks) 
presented in this paper extend beyond Tgff = 35,000 K. 

The synthetic colors have been tested and perform well in bandpasses equivalent to V or 
redder (with central wavelengths longer than ~ 5000 A) but the blue and ultraviolet bands 
(with central wavelengths shorter than ~ 5000 A) suffer from inaccuracy of the synthetic 
fluxes at shorter wavelengths. In cases where analysis in the bluer bands is important, 
empirical color transformations are strongly recommended. 



Cassisi et al.l (120041 ) explored the dependence of broadband synthetic colors on [a/Fe] for 



[a/Fe]=0 and +0.4. Their findings indicate that U-B and B-V are the colors most strongly 
effected by a-enhancement and that the influence of a-enhancement on broadband colors 
increases with increasing metallicity and decreasing Tes- The color transformations derived 
from the PHOENIX model atmosphere grid have a greater span of Tefr and [a/Fe], including 
[a/Fe] = -0.2 for the first time, and better resolution (0.2 dex) in [a/Fe]. The compar- 
isons perforrned wi th the PHOENIX synthetic colors qualitati vely agree with the findings of 
Cassisi et al.l (120041 ) though it is important to recognize that ICassisi et al.l (120041 ) compared 
colors at constant Z while in this paper comparisons are made at constant [Fe/H]. Hence 
an increase in [a/Fe] represents an increase in total Z and thus one expects the colors to 
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become redder with increasing [a/Fe] general rule. 

Figure [T] displays the changes to CMD morphology based solely on the level of a- 
enhancement at [Fe/H]=0 using a 10 Gyr scaled-solar isochrone. The figure shows [a/Fe] = 
0, +0.4, and +0.8 colors. The comparisons were assembled using a scaled-solar isochrone in 
order to demonstrate the effect of a-enhancement on the colors only. Two colors, V-I and 
J-Kg are shown in Figure [H The figure shows that V-I increases with [a/Fe] and while the 
difference is small near the MSTO, it grows as Tes decreases, becoming largest near My = 
10. The situation in J-Kg is reversed: the highest level of [a/Fe] produces the bluest J-Kg 
color while the lowest level produces the reddest color. The J-Kg result is, at first glance, 
counter-intuitive. However, the J band is particularly sensitive to absorption from water 
and, in the presence of enhanced oxygen (from [a/Fe]), substantially more water is able to 
form. This pushes fiux away from the J band and causes the J-Kg color to become bluer 
with increasing [a/Fe]. 



As reported by ICassisi et al.l (120041 ). the color differences due to variations in [a/Fe] 
decrease as total metallicity decreases. However, isochrones presented in this paper suggest 
that the separation along the lower MS in J-Kg remains noticeable even at [Fe/H]=-2 and 
this feature may prove a useful diagnostic for near infrared studies of metal poor systems of 
sufficient depth to see the lower MS. 

The remainder of this section give a brief description of the photometric systems into 
which the stellar evolution tracks and isochrones have been transformed. 



4.1. Johnson— Cousins UBVRI and the Two Micron All Sky Survey 



Besselll (Il990l ) defined filter transmission curves for the Johnson-Cousins UBV(RI)c 
photometric system . These filter curves were combined with 2MASS JHKg filters from 



Cohen et al.l (120031 ) to create a set of isochrones suitable for use with ground-based pho- 



tometry from the ultraviolet to the ne ar infrared. The UB VRI magnitudes were calibrated 
to Vega by following Appendix A of iBessell et al.l (119981 ) . The 2 MASS magnitudes were 
calibrated by setting the magnitudes of Vega to zero in each filter (ICohen et al.ll2003l ). 



4.2. The Hubble Space Telescope 



Bedin et al.l (120051 ) defined a method for generating synthetic photometry based on the 
detailed filter transmission curves for the Hubble Space Telescope ('HS T) Advanced Camer a 
for Surveys (ACS) and Wide Field Planetary Camera 2 (WFPC2) from lSirianni et al.l (120051 ). 
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Fig. 1. — The influence of ct-enhancement on broadband V-I and J-Kg colors at [Fe/H] = 
0. The underlying isochrones shown are all scaled-solar composition with [Fe/H] = at 10 
Gyr. Only the color transformations include a-enhancement. See the text for a discussion 
of the trends shown in the figure. 
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The synthetic photometry presented here are all on the Vegamag system. The following filters 
are included for the ACS-WFC: F435W, F475W, F555W, F606W, F625W, F775W, F814W, 
and F850LP. The following filters are included for the WFPC2: F336W, F439W, F555W, 
F606W, F791W, F814W, and F850LP. The performance of the sy nthetic color tran s forma - 



Dotter et al. 


(2007 


); 


Richer et al. 


(2007 


); 


Sieeel et al. 


(2007) 



4.3. The Sloan Digital Sky Survey 

SDSS produced (and its ongoing extensions continue to produce) vast quantities of stel- 
lar photometry and spectroscopy with important implications for our knowledge of galactic 
structure and evolution. The SDSS ugriz photometric system is based on the AB magnitude 
system but each of the five filters has an additional offset from AB of at most a few hun- 
dredths of a magnitude. The offsets, in the sense of mag(SDSS)-mag(AB) are: 0.04 in u, 
0.01 in g, in r, -0.015 in i, and -0.03 in z. The filter transmission curves and procedure 
for transforming synthetic fluxes to SDSS magnitudes magnitudes is outlined on the SDSS 
Photometric Flux Calibration web page [^. 

Typically, the synthetic colors provide good fits to SDSS data in g, r, and i. Synthetic 
photometry in the u band is subject to the same shortcomings as other bands at bluer 
wavelengths while problems with the underlying calibration of the SDSS photometric systems 
persist in both u and z. 



4.4. Semi— Empirical Transformations 

For general use, and comparison with the synthetic color transformations, t he tracks and 



isochrones are presented in colors based on the semi-empirical transformations of lVandenberg &: Clem 



(120031 . hereafter VC03) in BV(RI)c- VC03 compiled a set of Tefr-color transformations and 



„^ ^ 

(Bell & Gustafsson 


1978 


1989: ^ 


Vandenbere & Bell 


models (consistent with 


Castell 


1 (1999 


)) for higher 



temperatures. The color transformations are semi-empirical in the sense that VC03 began 
with a collection of synthetic colors and then applied corrections to the colors, primarily for 
cooler MS stars approaching and above [Fe/H]=0, based on observational constraints from 



http: / /www. sdss.org/dr6 /algorithms / fluxcal.html 
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open clusters, globular clusters, and field stars. 



Clem et al. (200 



system- 



Clem et al. 



) pre sented color transformations for the Stromgren uvby photometric 



(120041) used a similar approach to VC03 by adopting M ARCS models 



( Houdashelt et al. 200ol ) for Tefr < 8,000 K and colors based on Ku rucz models ( Castelli et al. 
1997h for higher temperatures. The uvby color transformations of I Clem et al.l (120041 ) provide 



much better fits to observations than any synthetic color tables (see that paper for a detailed 
discussion) and thus no synthetic Stromgren photometry has been provided. 



5. Comparisons 



This section will cor npare the DSEP mo dels to comparable models and observational 
data from the literature. iDotter et al.l (120071 ) presented comparisons with old, metal-poor 
isochrones (12 Gyr, Z=0.0001) from several groups and compared the DSEP isochrones 
to ACS photometry of t wo globular clusters: M 92 and 47 Tuc from the ACS Galactic 
Globular Cluster Survey (ISarajedini et al.l 120071 ). As a result, this section will focus on 
model comparisons at [Fe/H] > 0, along with younger, more metal rich open clusters. The 
one exception is a comparison of SDSS photometry from M 13 and NGC 2420 to the ugriz 
isochrones in S5.2.6. 



5.1. Comparisons with Other Models 



F ollowi ng Potter et al. ( 2007 ). comparison isoch rones were dra.wn fro m lPietrinferni et al 



(l2004j.l2006h:ICordier et all (120071 . hereafter BaSTII. ICirardi et all (120001 hereaf t er Pa dova), 
Vandenberg et al.l (120061 . hereafter Victoria- Regina (VR)), and lYi et al.l (120031 . 12004 here- 
after Yale-Yonsei (Y^)). In this section, the comparison is extended to solar metallicity and 
above as well as to younger ages. Comparisons are shown in both the log Tefj-log L/Lq plane 
and the V vs. V-I color-magnitude diagram. In each case, the DSEP isochrone is shown 
with the synthetic color transformation. A comparison of different color transformations 
applied to the DSEP isochrones appears later in this section. 

Figures [2] and [3] show isochrones at [Fe/H] ~+0.15 and an age of 625-630 Myr (com- 
patible with the the Hyades, see §5.2.2). The isochrones range in initial Z=0.025 to 0.03 
and Y=0.28 to 0.3. In Figure [2|, the Padova, BaSTI, and DSEP isochrones all show the 
AGE though the point at which the AGB terminates differs: DSEP and BaSTI stop at or 
before the onset of thermal pul sations though the BaSTI library has been updated to in- 
clude extended AGB evolution (ICordier et al.l 120071 ). Fainter than log L/L0=-1, the large 
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Fig. 2. — A plot of isochrones at [Fe/H] ~+0.15 and an age of ~ 625 Myr. The Padova, 
BaSTI, and DSEP isochrones show the He burning phase but the others do not. Good 
agreement exists among all the isochrones except for Padova that is cooler and fainter around 
the main sequence turn off. 
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Table 1. Minimum Mass (M0) for Convective Core Overshoot 
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1 KA'V 
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U.oo 
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[Fe/H] 


[a/Fe]=-0.2 
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1.6 
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1.5 


1.3 
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1.5 
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1.5 
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1.4 
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1.4 
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1.2 


1.2 
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1.0 


1.2 


1.2 


1.2 


1.2 


0.0 


1.2 


1.2 


1.1 


1.1 
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0.9 


1.2 


1.1 


1.2 


1.1 


+0.15 




1.1 


















+0.3 


1.1 


1.1 


1.1 
















+0.5 


1.1 


1.1 


1.0 


















Fig. 3. — The same as Figure [2] but focused on the region near the main sequence turn off. 
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Table 2. Mass Fractions of Z 



Element 




-0.2 













0.2 


[a/Fe]: 


0.4 






0.6 






0.8 




C 


2, 


.296E- 


■01 


1 


.721E- 


■01 


1 


.233E- 


■01 


8, 


.490E- 


■02 


5, 


.685E- 


■02 


3 


.731E- 


■02 


N 


6 


.726E- 


■02 


5, 


.041E- 


■02 


3 


.611E- 


■02 


2, 


.487E- 


■02 


1, 


.665E- 


■02 


1 


.093E- 


■02 





3, 


.940E- 


■01 


4, 


.680E- 


■01 


5, 


.314E- 


■01 


5, 


.800E- 


■01 


6, 


.155E- 


■01 


6 


.403E- 


■01 


Ne 


9, 


.044E- 


■02 


1, 


.050E- 


■01 


1, 


.165E- 


-01 


1, 


.271E- 


■01 


1, 


.349E- 


■01 


1 


.403E- 


■01 


Na 


2, 


.773E- 


■03 


2, 


.078E- 


■03 


1, 


.489E- 


■03 


1, 


.025E- 


■03 


6, 


.866E- 


■04 


4 


.506E- 


■04 


Mg 


3, 


.366E- 


■02 


3, 


.998E- 


■02 


4, 


.540E- 


■02 


4, 


.955E- 


■02 


5, 


.258E- 


■02 


5 


.470E- 


■02 


Al 


4, 


.814E- 


■03 


3 


.608E- 


■03 


2, 


.585E- 


■03 


1 


.780E- 


■03 


1 


.192E- 


■03 


7, 


.822E- 


■04 


Si 


3 


.715E- 


■02 


4, 


.412E- 


■02 


5, 


.OlOE- 


■02 


5, 


.468E- 


■02 


5, 


.803E- 


■02 


6 


.036E- 


■02 


P 


6, 


.493E- 


■04 


4, 


.866E- 


■04 


3 


.486E- 


■04 


2, 


.401E- 


■04 


1 


.607E- 


■04 


1 


.055E- 


■04 


S 


1 


.851E- 


■02 


2, 


.199E- 


■02 


2, 


.497E- 


■02 


2, 


.725E- 


■02 


2, 


.892E- 


■02 


3 


.008E- 


■02 


CI 


3 


.900E- 


■04 


2, 


.923E- 


■04 


2, 


.094E- 


■04 


1, 


.442E- 


■04 


9 


.656E- 


■05 


6 


.337E- 


■05 


Ar 


5, 


.793E- 


■03 


4, 


.342E- 


■03 


3, 


.lllE- 


■03 


2, 


.142E- 


■03 


1, 


.434E- 


■03 


9 


.414E- 


■04 


K 


3 


.045E- 


■04 


2, 


.282E- 


■04 


1, 


.635E- 


■04 


1, 


.126E- 


■04 


7, 


.539E- 


■05 


4 


.948E- 


■05 


Ca 


3 


.268E- 


■03 


3 


.882E- 


■03 


4, 


.408E- 


■03 


4, 


.811E- 


■03 


5, 


.106E- 


■03 


5 


.311E- 


■03 


Ti 


1, 


.519E- 


■04 


L 


.804E- 


■04 


2, 


.048E- 


■04 


2, 


.236E- 


■04 


2, 


.372E- 


■04 


2, 


.468E- 


■04 


Cr 


1 


.470E- 


■03 


1 


.102E- 


■03 


7, 


.894E- 


■04 


5, 


.436E- 


■04 


3 


.640E- 


■04 


2 


.389E- 


■04 


Mn 


1 


.075E- 


■03 


8, 


.054E- 


■04 


5, 


.771E- 


■04 


3 


.974E- 


■04 


2, 


.661E- 


■04 


1 


.746E- 


■04 


Fe 


1 


.020E- 


■01 


7, 


.641E- 


■02 


5, 


.474E- 


■02 


3 


.770E- 


■02 


2, 


.524E- 


■02 


1 


.657E- 


■02 


Ni 


6, 


.026E- 


■03 


4, 


.516E- 


■03 


3 


.235E- 


■03 


2, 


.228E- 


■03 


1 


.492E- 


■03 


9 


.791E- 


■04 



Table 3. Initial Compositions 







Y=0.245+1.54Z 






Y=0.33 


Y=0.4 


[Fe/H] [Q/Fe]=-0.2 





0.2 0.4 


0.6 


0.8 


0.4 


0.4 



-2.5 


"V 

A= 


=0.7549 





TC A O 

7548 





7548 


0.7547 


0.7546 


0.7543 


0.DD99 


0.6699 


0.6000 


A C AAA 

0.5999 




z= 


4.09E-5 


5 


48E-5 


7 


63E-5 


l.llE-4 


1.65E-4 


2.52E-4 


4.85E-5 


9.85E-5 


4.34E-5 


8.82E-5 


-2.0 


x= 


=0.7547 





7545 





7544 


0.7541 


0.7536 


0.7529 


0.6698 


0.6697 


0.5999 


0.5997 




Z= 


1.30Iii-4 


1 


72E-4 


2 


41E-4 


3.50E-4 


5.24E-4 


7.97E-4 


1.53E-4 


O C AT71 A 

3.50E-4 


1.37E-4 


2.79E-4 


-1.5 


X- 


=0.7539 





7536 





on 

7530 


0.75/1 


0.7507 


0.7484 


0.DD95 


0.6690 


A C A A/^ 

0.5996 


A A A 1 

0.5991 




z= 


4.10E-4 


5 


47E-4 


7 


62E-4 


l.llE-3 


1.65E-3 


2.50E-3 


4.85E-4 


9.83E-4 


4.34E-4 


8.81E-4 


-1.0 


x= 


=0.7516 





7505 





7487 


0.7459 


0.7415 


0.7346 


0.6685 


0.6669 


0.5986 


0.5972 




z= 


1.29E-3 


1 


72E-3 


2 


40E-3 


3.46E-3 


5.15E-3 


7.77E-3 


1.53E-3 


3.46E-3 


1.37E-3 


2.78E-3 


-0.5 


x= 


=0.7444 





7441 





7355 


0.7270 


0.7139 


0.6942 


0.6652 


0.6603 


0.5957 


0.5913 




z= 


4.05E-3 


5 


37E-3 


7 


44E-3 


1.07E-2 


1.57E-2 


2.32E-2 


4.82E-3 


9.70E-3 


4.31E-3 


8.69E-3 


0.0 


x= 


=0.7174 





7071 





6882 


0.6628 


0.6247 


0.5710 


0.6550 


0.6402 


0.5866 


0.5734 




z= 


1.43E-2 


1 


89E-2 


2 


55E-2 


3.52E-2 


4.97E-2 


7.02E-2 


1.50E-2 


2.98E-2 


1.34E-2 


2.66E-2 


+0.15 


x= 







6880 






















z= 




2 


56E-2 




















+0.3 


x= 


=0.6842 





6638 





6334 


















z= 


2.70E-2 


3 


48E-2 


4 


64E-2 
















+0.5 


x= 


=0.6491 





6185 





5792 


















z= 


4.04E-2 


5 


21E-2 


6 


71E-2 

















Table 4. Surface Compositions for ZAHB Models 









Y™t=0.2454 


-1 54Z 










=0.33 


V- ■ 


=0.4 


fFe/Hl 


[a/Fe]=-0.2 







0.2 


0.4 


0.6 


0.8 







0.4 





0.4 


-2.5 


X=0.7564 





7565 


0.7561 


0.7549 


0.7549 


0.7546 





6835 


0.6828 


0.5967 


0.6212 




Z=3.97E-5 


5 


31E-5 


7.40E-4 


1.60E-4 


1.08E-4 


2.44E-4 


4 


63E-5 


9.40E-5 


8.40E-3 


8.33E-5 


-2.0 


X=0.7547 





7548 


0.7551 


0.7541 


0.7524 


0.7510 





6820 


0.6810 


0.6208 


0.6193 




Z=1.26E-4 


1 


67E-4 


2.33E-4 


3.41E-4 


5.06E-4 


7.69E-4 


1 


46E-4 


2.97E-4 


1.30E-4 


2.64E-4 


-1.5 


X=0.7529 





7517 


0.7509 


0.7491 


0.7464 


0.7427 





6785 


0.6763 


0.6180 


0.6156 




Z=3.96E-4 


5 


29E-4 


7.37E-4 


1.07E-3 


1.60E-3 


2.43E-3 


4 


64E-4 


9.44E-4 


4.12E-4 


8.37E-4 


-1.0 


X=0.7472 





7453 


0.7423 


0.7383 


0.7325 


0.7238 





6732 


0.6698 


0.6125 


0.6076 




Z=1.25E-3 


1 


67E-3 


2 33E-3 


3.37E-3 


5.02E-3 


7.59E-3 


1 


47E-3 


3.00E-3 


1.31E-3 


2.67E-3 


-0.5 


X=0.7350 





7309 


0.7242 


0.71485 


0.7010 


0.6810 





6641 


0.6600 


0.6034 


0.5967 




Z=3.94E-3 


5 


24E-3 


7.27E-3 


1.05E-2 


1.54E-3 


2.28E-2 


4 


67E-3 


9.48E-3 


4.16E-3 


8.41E-3 


0.0 


X=0.7166 





6928 


0.6741 


0.6497 


0.6131 


0.5640 





6477 


0.6303 


0.5877 


0.5854 




Z=1.41E-2 


1 


85E-2 


2.50E-2 


3.46E-2 


4.88E-2 


6.87E-2 


1 


47E-2 


2.91E-2 


1.31E-2 


2.68E-2 


+0.15 


x= ... 
z= ... 



2 


6737 
51E-2 




















+0.3 


X=0.6696 
Z=2.66E-2 



3 


6505 
42E-2 


0.6204 
4.57E-2 


















+0.5 


X=0.6360 
Z=3.98E-2 




5 


6075 
12E-2 


0.5734 
6.57E-2 



















^Yi„it and Zi„it refer to the initial Y and Z values, for the initial compositions see Table |3l 
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differences are due to different mass ranges, the EOS, and low temperature opacities used 
by each group. The BaSTI isochrone begins at 0.5 Mq, VR and isochrones begin at 0.4 
Mq, the Padova isochrone begins at 0.15 Mq, and the DSEP isochrone begins at 0.3 Mq at 
this age. 

In Figure [3] the agreement is quite good among DSEP, BaSTI, and isochrones. The 
Padova isochrone is less luminous in the region of the MS and subgiant branch (SGB) while 
the VR isochrons is more luminous in the same region. These differences in luminosity are 
primarily due to different adopted treatment of convective core overshoot. 

Figure [2] indicates that the Y^ isochrone extends to cooler temperatures than any of 
the others despite the fact that it terminates at 0.4 Mq while DSEP and Padova terminate 
at considerably lower mass. This effect is almost certainly due to the EOS used in the Y^ 
calculations. 

Figures m and [5] show isochrones at solar metallicity and 4 Gyr (compatible with M 67, 
see §5.2.4). Each isochrone has an initial Z=0.018-0.020 and Y=0.27-0.28; the adopted 
physics (especially the EOS and convective core overshoot prescription) and solar heavy 
element mixture differ from one group to another. In the left panels of Figure H] and O 
the isochrones are in reasonable agreement for points brighter than log L/Lq=-1 (Padova 
and BaSTI isochrones include He burning stars while the others do not). In Figure O the 
BaSTI, VR, and DSEP isochrones all agree well on the temperature and luminosity of the 
MS turn off (TO) and SGB. In the right panels of Figure H] and 0, all isochrones agree well at 
magnitudes brighter than My=6 except for Padova that is consistently redder and brighter 
on the upper RGB. In Figure H] the DSEP isochrone includes the complete lower mass range 
(down to 0.1 Mq) and the fact that the Y^ isochrone is cooler and almost as faint in My is 
significant. Conversely, the Padova isochrone is hotter and bluer than the others for points 
fainter that log L/Lq=-1 (My=7) although the mass-luminosity relation is similar to the 
DSEP isochrone. 

Figures El [3], H] and [5] show general agreement amongst the different sets of isochrones. 
The main exceptions are the Padova isochrones that appear cooler and redder near the 
MSTO and on the RGB but hotter and bluer on the lower main sequence. Amongst the 
other models the largest differences occur near the MSTO where the adopted core overshoot 
treatments differ and on the lower main sequence where the adopted EOS and minimum 
masses differ. The DSEP models use the same EOS and opacities as the VR and BaSTI 
models and it is reassuring that these models all lie close together aside from the noted 
exceptions due to the included mass range and treatment of convective core overshoot. 
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Fig. 4. — The figure shows solar metalhcity isochrones at 4 Gyr. The Padova and BaSTI 
isochrones include the He burning phase but the others do not. There is remarkable agree- 
ment in the H-R diagram between the tip of the red giant branch and log L/Lq=-1. Below 
this point, the range of masses plotted and the adopted EOS cause large differences. There 
is good agreement in the color-magnitude diagram between Mv=6 and but beyond this 
the mass range plotted, the underlying physics, and the choice of color transformation cause 
significant discrepancies. 
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Table 5. He Core Mass (Mq) for ZAHB Models 









=0.245+ 


1.54Z 






Y= 


0.33 




0.4 


fFe/Hl 

L / J 


[a/Fe]=-0.2 





0.2 


0.4 


0.6 


0.8 





0.4 





0.4 


-2.5 


0.506 


0.505 


0.504 


0.502 


0.500 


0.497 


0.489 


0.485 


0.474 


0.473 


-2.0 


0.500 


0.499 


0.498 


0.494 


0.492 


0.491 


0.483 


0.479 


0.469 


0.468 


-1.5 


0.495 


0.493 


0.491 


0.489 


0.488 


0.486 


0.478 


0.474 


0.464 


0.463 


-1.0 


0.489 


0.488 


0.486 


0.484 


0.481 


0.478 


0.473 


0.469 


0.460 


0.458 


-0.5 


0.483 


0.482 


0.480 


0.477 


0.473 


0.469 


0.469 


0.465 


0.457 


0.455 


0.0 


0.475 


0.473 


0.469 


0.464 


0.458 


0.447 


0.465 


0.460 


0.454 


0.451 


+0.15 




0.468 


















+0.3 


0.468 


0.464 


0.459 
















+0.5 


0.461 


0.456 


0.449 


















Fig. 5. — The same as Figure H] but focused on the region near the main sequence turn off. 
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5.2. Comparisons with Observational Data 



Several papers have demonstrated the performance of the DSEP isochrones when apphed 
to HST/ACS data in two of the most popu l ar bands for globular cluste r photometry ( F606W 



and F814W) including: ISaraiedini et all (120071 ). iDotter et al.l (120071 ). ISiegel et al.l (120071 ) 



and [Richer et al.l (120071 ). This section will focus on the ground-based photometric systems 
presented in this paper and also compares the synthetic color transformation with the VC03 
transformation in B-V and V-I. 



Grocholski fc Sarajedinil (120031 . hereafter GS03) combined optical photometry of six 
Galactic open clusters with near- infrared data fro m 2MASS. Similar compilations for the 
Hyades (jPinsonneault et al.l 12004 ) and NGC 6791 (ICarney et al.ll2005l ) are presented along 
with three clusters from GS03 to span a range of ages and metallicities in open clusters. 
These combined optical and near-infrared data provide an excellent source of comparisons 
with the DSEP isochrones. Table E] lists the relevant parameters of the clusters featured in 
this section. 



The adopted reddening values are listed as E(B-V). The extinction curve of I Johnson 

(119651 ) was used to determine the reddening in other filter combinations: E(V-I) = 1.4 



E(B-V), E(V-J) = 2.4 E(B-V) and E(V-Ks) = 2.8 E(B-V). 

These comparisons are not meant to determine the ages of the clusters but rather to 
demonstrate the performance of the isochrones. The ages were chosen to be consistent with 
other published results but if the age shown is considerably different from the true age of 
the cluster it will influence the quality of the fit. 



5.2.1. M 37 



Figure [6] shows the CMD of M 37 in B-V (IKalirai et al.ll200lh and V-J and V-Ks (GS03) 



along with a [Fe/H]=0, 500 Myr isochrone. The isochrone is plotted using the synthetic (solid 
line) and VC03 (dashed line) colors. 

In B-V (left panel of Figure [HD, the synthetic color isochrone is redder than the data 
by ~0.1 mag throughout the CMD while the VC03 isochrone provides a reasonable fit to 
the data from the MSTO to V=16 before drifting to the red of the data. The He burning 
portion of the isochrone is redder than the data by 0.1 mag in both color transformations. 

In V-J and V-Kg, the synthetic color isochrone provides an excellent fit to the data for 
the entire extent of the MS (the MS data cuts off at V~20 in V-J and V-Ks). The isochrone 
is redder than the data for the He burning stars in both V-J and V-Kg. 
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Fig. 6.— CMDs of M 37 from iKalirai et all (120011 ) and GS03 are plotted with a [Fe/H]=0, 
500 Myr isochrone. The isochrone is shown in synthetic (sohd hne) and VC03 (dashed hne) 
colors. The synthetic color isochrone is redder than the data in B-V but a good fit in V-J 
and V-Kg. The He burning portion of the isochrones is redder than the data in all cases. 
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We note in passing that the VC03 B-V isochrone shown in the left panel of figure [6] does 
not match the shape of the CMD for V > 17. Although M 37 and M 6 7 are both roughly 
[Fe/H ] = clusters their B-V CMDs from lKalirai et al.l J200l[ ) and M 67 jMontgomerv et al. 



19941 ). respectively, adjusted for distance and reddening, appear to disagree somewhat for 
My > 6 when viewed together. Since the VC03 colors are based in part on the CMD of 
M 67, it is impossible for them to fit both clusters simultaneously. 



5.2.2. The Hyades 

Figure [7| shows the Hipparcos B-V (left panel) and V-I (center) CMDs of the Hyades 



(Ide Bruiine et al.ll200l[) a nd the combined Hipparcos and 2MASS V-Kg CMD compiled by 



Pinsonneault et al.l (120041 ). Along with the data, a [Fe/H] =+0.15, 625 Myr isochrone is 
shown in the synthetic (solid line) and VC03 (dashed line) color transformations. The 
choice of a 625 Myr isochrone may be an underestimate of the age which is typically quoted 
as between 650 and 700 Myr (e.g. VC03). The He abundance assumed for the [Fe/H]=+0.15 
models may also be a cause for disagreement since the adopted AY/AZ yields Yj„jf =0.2864 
whereas the Hyades Y value inferred from model comparisons is lower by ~0.02 (see VC03 
and references therein for a discussion). 

In B-V, the synthetic color isochrone is redder than the data by 0.05-0.1 mag down to 
My=7 where the slope changes and the isochrone becomes bluer than the data by Mv=8. 
The VC03 color isochrone, on the other hand, is a good fit for the entirety of the MS. Both 
isochrones appear cooler than the two He burning stars in the B-V CMD. 

In V-I, the color transformations are quite similar and both perform well down to My=6. 
At this point, the isochrones become redder than the data before matching the data again 
at My=8. 

In V-Ks, the synthetic color isochrone fits the data for stars with My > 3 but becomes 
bluer than the data just below the MSTO. The discrepancy just below the MSTO would be 
avoided by plotting a slightly older isochrone. 



5.2.3. NGC 2420 

Figure E shows the B-V, V-I, and V-Kg CMDs of NGC 2420 from GS03. Plotted 
along with the data is an isochrone interpolated to [Fe/H]=-0.25 and 2.15 Gyr. Given the 
adopted distance and reddening in Table [6], the best fit age plotted in Figure [8] is older than 
the age quoted by GS03. This is not surprising since ultimately the GS03 age estimate for 




Fig. 7. — The Hyades in B-V and V-I from Ide Bruijne et al.l (l200lh along with V-Kg from 
Pinsonneault et al.l (120041 ) . Along with the data, an isochrone with [Fe/H]=+0.15 at 625 
Myr is plotted using the synthetic colors (solid line) and VC03 colors (dashed line). 
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Fig. 8. — NGC 2420 from GS03 along with an isochrone interpolated to [Fe/H]=-0.25 and 
2.15 Gyr. The isochrone is shown with synthetic colors (solid line) and VC03 colors (dashed 
line). Also shown is a 1.5 Mq He burning track interpolated to [Fe/H]=-0.25. 
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NGC 2420 was based on Padova isochrones and these were shown to be fainter and redder 
near the MSTO than the DSEP isochrones (see Figures [2] and H] and the discussion in §5.1). 
As a result, the best fit Padova isochrone will be younger than the best fit DSEP isochrone. 

In B-V the synthetic isochrone is redder than the data, particularly on the lower MS 
and the RGB, though it matches well near the MSTO. The VC03 color isochrone fits the 
entire CMD well. The He burning track lies fainter than the red clump for both color 
transformations with the color predicted correctly by VC03 but too red by the synthetic 
color. If NGC 2420 is younger than 2.1 Gyr, then the stars in the red clump will be more 
massive and more luminous. 

In V-I, both synthetic and VC03 colors fit the CMD well except for the synthetic color 
on the RGB that is slightly cooler than the data. The color of the red clump is correctly 
predicted by VC03 but slightly too red in the synthetic color. 

In V-Ks, the synthetic isochrone fits the SGB and RGB and correctly predicts the color 
of the red clump. However, the V-Kg color appears to be too blue on the MS. 



5.2.4. M 67 



Figure 1 shows the B-V, V-I, and V-Ks CMDs of M 67 from Montgomery et all fll994f l 
and GS03 along with a [Fe/H]=0, 4 Gyr isochrone. The isochrone has been transformed 
using the synthetic (solid line) and VC03 (dashed line) colors. A He burning track with 
[Fe/H]=0 and M=1.25 Mq is shown to indicate the predicted location of the red clump. 

The fit of the synthetic color isochrone to the data in B-V (left panel of Figure [9]) is 
poor, with the isochrone redder than the data by about 0.05-0.1 mag throughout the CMD 
while the VC03 color isochrone provides an excellent fit. The He burning model has the 
correct V magnitude but is slightly redder than the observed red clump in B-V for both 
color transformations. 

In V-I (middle panel of Figure [9]), both synthetic and VC03 colors match the data well 
over most of the CMD. The only exception is that, for stars fainter than V=17, the VC03 
colors are redder than the data while the synthetic colors follow the data closely down to 
V=20. The He burning model correctly predicts the color of the red clump in V-I for both 
color transformations. 

In V-Kg only the synthetic color isochrone is plotted in the right panel of Figure [91 
As with V-I, the synthetic V-Kg color is an excellent fit to the data for the lower MS, the 
MSTO, the RGB, and the location of the red clump. 
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Fig. 9.— The CMDs of M 67 in B-V and V-I flMontgomerv et al.lll994l ) and V-Ks (GS03). 
Isochrones with [Fe/H]=0 at 4 Gyr are shown using synthetic (sohd hne) and VC03 (dashed 
hne) colors. A He-burning track with M=1.25 Mq at [Fe/H]=0 is plotted to indicate the 
predicted level of the red clump. 
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5.2.5. NGC 6791 



Carney et al.l (I20051) combined JHK photometry of NGC 6791 with optical data from 
Stetson et al.l (120031 ) . The resulting CMDs in B-V, V-I, and V-K are shown in Figure [TOl 
The V-K CMP has been adjusted by +0.024 mag to convert from CIT K to 2MASS 
(ICarpenterll200ll ). The isochrones shown in Figure [TU] have [Fe/H].^^^ = +0.36 and an age of 
9 Gyr. 

The B-V CMD continues to be problematic for the synthetic colors where, using the 
adopted E(B-V), the synthetic isochrone is ~0.1 mag redder than the data. The isochrone 
with VC03 colors is a substantially better fit to the da ta. This is no t surprising since the 
VC03 colors were calibrated at high metallicity with the lStetson et al.l (120031 ) data. The V-I 
CMD is equally well fit by both synthetic and VC03 colors. The synthetic color isochrone 
in V-K is an excellent fit to the shape of the CMD as well. 

For the isochrone used in Figure [101 neglecting mass loss, the mass at the tip of the 
RGB is 1.12 Mq. a 1 Mq He burning track is plotted to show the predicted level of the 
red clump in Figure [TOl The He burning track provides a good fit to the V mag of the red 
clump but lies a few hundredths of a magnitude redder than the red clump in B-V for both 
color transformations. In V-I, the color of the red clump is in excellent agreement with the 
model. In V-K, the model is bluer than the red clump by a few hundredths of a magnitude. 



It should be noted that while the adopted reddeni ng value is EfB-V) = O.l JCarney et al. 



(|2005[ ) found E(B-V) = 0.14 at [Fe/H]=+0.4 and the ISchlegel et al.l (Il998h dust maps yield 
E(B-V) = 0.155 for NGC 6791. These values are both higher than the adopted value of 
E(B-V) and would alter the isochrone fits shown in Figure [TDl 

If instead E(B-V) = 0.15 is adopted as the reddening then the isochrone fits re- 
quire either younger age or lower metallcity (or both). Using the same isochrones (with 
[Fe/H]^^^^=+0.36) but with the larger reddening, a reasonably good fit is achieved at 8 Gyr 
but with an increased distance modulus of DMy = 13.45. Reducing [Fe/H] to +0.2 intro- 
duces a substantial mismatch between the observed CMD and the models along the MS with 
the isochrones becoming bluer than the data at 2 mag below the MSTO. 

Given that the adopted [Fe/H] value provides a good fit to the entir e CM D and agrees 



with recent spectroscopic estimates by, for example, ICarraro et al.l (120061 ) and lOriglia et al. 



(120061 ). the isochrones show a clear preference for the lower reddening value of E(B-V)=0.1 
used in Figure [TOl Constrainin g the distance, reddeni ng, composition, and age of NGC 6791 
remains a thorny problem (see IChaboyer et al.l (119991 ) for a detailed discussion). Because of 
it high metallicity, a conclusive age determination of NGC 6791 will most likely require that 
the models follow the detailed abundance patterns found by spectroscopy and not simply 
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Fig. 10.— The B -V, V-I, V-K CMDs of NGC 6791 from IStetson et all (j2003|) and 
Carney et al.l (120051 ). Plotted on top of the data are isochrones with [Fe/H]=+0.36 at 9 
Gyr in synthetic (sohd line) and VC03 (dashed hne) colors. A 1 M© He burning track of the 
same composition is plotted to indicate the predicted location of the red clump. 
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match the observed [Fe/H]. For example, lOrigha et aLl (120061 ) found that [a/Fe] 
[C/Fe]=-0.35 in NGC 6791. 



but 



5.2.6. Clusters in the SDSS ugriz system 



Lee et al.l (120071 ) present spectroscopically selected cluster members for stars in a sample 



of galactic globular and open clusters for purposes of validating the SEGUE Stellar Parameter 
Pipeline. The selected stars from M 13 and NGC 2420 are used here because they represent 
the cleanest available data sets of clusters on the SDSS ugriz system. Though these sequences 
are quite narrow and well-defined, they are not as deep as allowed by SDSS, but are suitable 
for the purpose of demonstrating the performance of the isochrones. 

Figure [n] shows the g vs. (g-i)o CMDs of M 13 (left) and N GC 2420 (r i ght) a long with 
appropriate isochron es. The photometry has been de- reddened by lLee et al.l (120071 ) using the 
Schlegel et al.l (1l998l ) dust maps. Plotted on top of the M 13 data in Figure fTTl are isochrones 
with [Fe/H]=-1 .5, [Q;/Fe]=+ 0.2 at 11 and 13 Gyr assuming a distance modulus of 14.4 in 
agreement with iHarrij (119961 ). The isochrones follow the locus of stars on the MS and RGB 
well. The only real discrepancy is the slope of the SGB but this is a common problem with 
synthetic color isochrones compared to metal poor systems. 

Plotted on top of the NGC 2420 data in Figure ITT] is the same isochrone as in §5.2.3 but 
now transformed to g vs. g-i plane. The agreement between the data and the isochrone is 
excellent on the MS but the observational constraints of the sample do not permit comparison 
with stars in the convective hook or past the MSTO. 



5.3. Discussion 

Bluer synthetic colors, as represented by B-V in this section, continue to require adjust- 
ments in order to match observed CMDs as Figures [6]through[T0l indicate. The semi-empirical 
color transformation of VC03 are an improvement over the synthetic B-V color but seem 
unnecessary for V-I when used in conjuction with the DSEP isochrones. While the synthetic 
colors and magnitudes are suitable analysis of observational data in bandpasses with central 
wavelengths longer than ~ 5000 A, empirically adjusted colors ought to be used at shorter 
wavelengths. In these redder bands, the synthetic colors perform well throughout the CMD, 
including on the lower main sequence where they perform better than the empirical colors. 

The optical to near-infrared comparisons in V-J (M 37 only) and V-K of the open 
clusters are favorable for entire extent of the MS and the RGB in the older clusters but 
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Fig. 11.— The de-reddened g-i CMDs of M 13 and NGC 2420 from lLee et all (120071 ) plotted 
along with isochrones of appropriate age and metallicity transformed to the SDSS ugriz 
system. 
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larger reddening than suggested by the adopted reddening curve is needed in some cases. 

The location of He burning stars in the open clusters is matched by the models in the 
older clusters (M 67 and NGC 6791) but for the younger clusters the models predict the 
red clump to be fainter (NGC 2420) and/or redder (M 37 and the Hyades) than the data 
suggest. The general lack of evolved stars in the CMDs of both M 37 and the Hyades makes 
it difficult to speculate as to whether the models are intrinsically too cool or whether the 
colors are too red for a given Teg. 



6. Summary 

The Dartmouth Stellar Evolution Program has been used to compute a larg e grid of stel 



l ar evo lution tracks. A subset of these tracks and isochrones were presented by iDotter et al. 



( 120071 ) and are extended in this paper to include younger ages, higher metallicities, and 
transformations to multiple photometric systems. These models take advantage of the latest 
advances in radiative opacities, nuclear reactions rate, and the equation of state. 

The Dartmouth Stellar Evolution Database comprises a set of stellar evolution tracks 
and isochrones transformed to several photometric systems along with a set of computer 
programs that allow for interpolation and construction of luminosity functions and synthetic 
horizontal branch models. The database spans a wide range of compositions from [Fe/H]=- 
2.5 to +0.5, [a/Fe]=-0.2 to +0.8 (at [Fe/H] <0) or +0.2 (at [Fe/H] >0), and additional 
models with enhanced He for [Fe/H] =-2.5 to 0. 

The models were transformed to the observational plane using a sy nthetic color trans- 



forma tion based on PHOENIX model atmospheres for Teg <10,000 K and ICastelli fc Kurucz 



( 120031 ) for hotter temperatures. The synthetic color transformations are compared to the 
VC03 semi-empirical transformations in B-V and V-I. The synthetic colors perform well 
in V-1 but are consistently too red in B-V. As a general rule, the synthetic colors should 
only be trusted in bandpasses equivalent to V or redder, i.e. those with central wavelengths 
longer than ~ 5000 A. In these bands, the synthetic colors were found to perform well — and 
better than the empirical colors — on the lower main sequence. Synthetic photometry in the 
bluer bandpasses (those with central wavelengths shorter than ~ 5000 A) does not compare 
favorably with observations. Empirically adjusted magnitudes and colors should be favored 
for the bluer bandpasses whenever possible. 

The isochrones were compared to isochrones from other groups at [Fe/H]=0 and +0.15 
as well as photometry from open clusters and one globular cluster. The isochrones extend 
to lower mass (0.1 Mq) than those of other major isochrone libraries and provide reasonable 
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fits to the MS of the clusters shown in §5. 

The Dartmouth Stellar Evolution Database should prove useful in studies of resolved 
stellar populations in the local Universe and is suitable for application to population synthesis 
and integrated light models. The database is available through a dedicated website. 
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Table 6. Comparison Clusters 



ID 


BMv 


E(B-V) 


[Fe/H] 


Approximate Age (Gyr) 


Photometry 


Reference (s) 


M 37 


11.55 


0.27 


0.09 


0.5 


B V J Ks 


1,2 


Hyades 






0.13 


0.65 


B V I Ks 


3,4 


NGC 2420 


12.1 


0.05 


-0.27 


2 


B V I Ks 


1 


M 67 


9.65 


0.04 





4 


B V I Ks 


1,5,6 


NGC 6791 


13.35 


0.1 


0.4 


9 


B V I Ks 


7,8,9 



References. 



1-GS 03. 2-lKalirai et a. 



Pinsonneault et al.l (1200411. 5-lMontgomerv et al. 



Stetson et all J2003h !9 lciiabover et all Jl999h 



(1200 ih . 3- lde Bruiine et all pOOlh. 4- 
Jl994h . 6-tavlor. 7- fCarnev et all J2005h . 8- 



